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ABSTRACT 

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons partici¬ 
pating in large-scale bulk flows around central galaxies (CGs) at redshift z ~ 0.1. We estimate the pairwise momentum of the kSZ temperature 
fluctuations at the positions of the Central Galaxy Catalogue (CGC) samples extracted from Sloan Digital Sky Survey (SDSS-DR7) data. For the 
foreground-cleaned SEVEM, SMICA, NILC, and COMMANDER maps, we find 1.8-2.5 cr detections of the kSZ signal, which are consistent with the kSZ 
evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W band (3.3 cr). We further reconstruct 
the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctu¬ 
ations and estimates of CG radial peculiar velocities. This correlation function yields a 3.0-3.7 cr detection of the peculiar motion of extended gas 
on Mpc scales in flows correlated up to distances of 80-100 IP 1 Mpc. Both the pairwise momentum estimates and the kSZ temperature-velocity 
field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of > 1 Mpc, more 
than twice the mean virial radius of halos. This is consistent with the predictions from hydrodynamical simulations that most of the baryons are 
outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a 
semi-analytic model built upon A-body simulations, and interpret the proportionality constant as an effective optical depth to Thomson scattering. 
We find tt = (1.4 ± 0.5) x 1CT 4 ; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes 
little signal to X-ray or thermal Sunyaev-Zeldovich observations. 

Key words. Cosmology: observations - cosmic microwave background - large-scale structure of the Universe - Galaxies: clusters: general - 
Methods: data analysis 


1. Introduction 

The kinetic Sunyaev-Zeldovich effect (hereafter kSZ, Sunyaev 
& Zeldovich 1972, 1980) describes the Doppler boost expe¬ 
rienced by a small fraction of the photon bath of the cosmic 
microwave background (CMB) radiation when scattering off a 


cloud of moving electrons. In the limit of Thomson scattering, 
where there is no energy exchange between the electrons and the 
CMB photons, the kSZ effect is equally efficient for all frequen¬ 
cies and gives rise to relative brightness temperature fluctuations 
in the CMB that are frequency independent. The exact expres- 
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sion for this effect was first written as equation 15 of Sunyaev & 
Zeldovich (1970), and reads 

ST T (v \ 

y(n) = ~ J d/o-jn^- ■ ftj. ( 1 ) 

In this expression, the integral is performed along the line of 
sight, n e denotes the electron number density, cr x is the Thomson 
cross-section, and (v/c) ■ h represents the line of sight component 
of the electron peculiar velocity in units of the speed of light 
c. The equation above shows that the kSZ effect is sensitive to 
the peculiar momentum of the free electrons, since it is propor¬ 
tional to both their density and peculiar velocity. A few previous 
studies (Hernandez-Monteagudo et al. 2006b; Bhattacharya & 
Kosowsky 2008; Ma & Zhao 2014; Keisler & Schmidt 2013; Li 
et al. 2014) have proposed that it be used to trace the growth 
of velocities throughout cosmic history and its connection to 
dark energy and modified gravity. Zhang & Stebbins (2011) and 
Planck Collaboration Int. XIII (2014) have also used the kSZ 
effect to test the Copernican Principle and the homogeneity of 
the Universe. In addition, the impact of the kSZ effect on sub¬ 
cluster scales has also been investigated (Inogamov & Sunyaev 
2003; Dolag & Sunyaev 2013). 

With these motivations, there have been previous attempts to 
detect the kSZ effect in existing CMB data (Kashlinsky et al. 
2008, 2010; Lavaux et al. 2013; Hand et al. 2012). The re¬ 
sults of Kashlinsky et al. (2008, 2010), which point to the ex¬ 
istence of a bulk flow of large amplitude (800-1000 km s -1 ) ex¬ 
tending to scales of at least 800 Mpc, have been disputed by 
a considerable number of authors (e.g., Keisler 2009; Osborne 
et al. 2011; Mody & Hajian 2012; Planck Collaboration Int. 
XIII 2014; Feindt et al. 2013). On the other hand, Lavaux et al. 
(2013) have claimed the detection of the local bulk flow (within 
80/r 1 Mpc) by using the kSZ in WMAP data in the direction of 
nearby galaxies. While these results remain at a low (roughly 
1.7 cr) significance level, they are also in slight tension with 
the results that we present here. On the other hand, the work 
of Hand et al. (2012) constitutes the first clear detection of the 
kSZ effect (using the “pairwise” momentum approach that we 
describe below), and no other group has confirmed their re¬ 
sults to date. In addition, Sayers et al. (2013) have provided 
a first claimed detection of the kSZ effect in a single source. 
After the first data release of the Planck mission, some dis¬ 
agreement has been claimed between the cosmological frame 
set by Planck and measurements of peculiar velocities as in¬ 
ferred from redshift space distortions (Macaulay et al. 2013). 
More recently, Mueller et al. 2014 forecasted the detectability of 
the neutrino mass with the kSZ pairwise momentum estimator 
by using Atacama Cosmology Telescope (ACTPol) and Baryon 
Oscillation Spectroscopic Survey (BOSS) data. 

This work represents the second contribution of the Planck 1 
collaboration to the study of the kSZ effect, in which we focus 
on constraining the “missing baryons” with this signal (DeDeo 
et al. 2005; Hernandez-Monteagudo & Sunyaev 2008; Ho et al. 
2009; Hernandez-Monteagudo & Ho 2009; Shao et al. 2011). 
Numerical simulations (Cen & Ostriker 2006) have shown that 
most of the baryons lie outside galactic halos, and in a diffuse 


1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


phase, with temperature in range of 10 5 -10 7 K. This “warm-hot” 
intergalactic medium is hard to detect in X-ray observations due 
to the relatively low temperature. Thus, most of the baryons are 
“missing” in the sense that they are neither hot enough (T < 
10 8 K) to be seen in X-ray observations, nor cold enough (T > 
10 3 K) to be made into stars and galaxies. However, as we can 
see in Eq. (1), the kSZ signal is proportional to the gas density 
and peculiar velocity, and thus the gas temperature is irrelevant. 
Therefore the kSZ effect has been proposed as a “leptometer,” 
since it is sensitive to the ionized gas in the Universe while most 
of the baryons remain undetected. 

In this work, we measure the kSZ effect through two dis¬ 
tinct statistics: the kSZ pairwise momentum, which was used in 
Hand et al. (2012) to detect kSZ effect for the first time; and 
the cross-correlation between the kSZ temperature fluctuations 
with the reconstructed radial peculiar velocities inferred from a 
galaxy catalogue. These analyses allow us to probe the amount 
of gas generating the kSZ signal and provide direct evidence of 
the elusive missing baryons in the local Universe. In Sect. 2 we 
describe the CMB data and the Central Galaxy Catalogue that 
will be used in our analysis. In Sect. 3 we describe two statisti¬ 
cal tools we use. Section 4 presents and explains our results and 
illustrates the physical meaning of the significance. The conclu¬ 
sions and discussion is presented in the last section. Throughout 
this work we adopt the cosmological parameters consistent with 
Planck Collaboration XIII (2015): Q m = 0.309; O a = 0.691; 
n s = 0.9608; erg = 0.809; and h = 0.68, where the Hubble con¬ 
stant is Hq — 100/zkms -1 Mpc -1 . 

2. Data description 

2.1. Planck data 

This work uses Planck data that are available publicly, both 
raw frequency maps and the CMB foreground cleaned maps. 2 
The kSZ effect should give rise to frequency-independent tem¬ 
perature fluctuations and hence constitutes a secondary ef¬ 
fect with identical spectral behaviour as for the intrinsic CMB 
anisotropies. Therefore the kSZ effect should be present in all 
CMB frequency channels and in all CMB foreground-cleaned 
maps. However, one must take into account the different effec¬ 
tive angular resolutions of each band, particularly when search¬ 
ing for a typically small-scale signal such as the kSZ effect. The 
raw frequency maps used here are the LFI 70 GHz map, and the 
HFI 100, 143, and 217 GHz maps. These have effective FWHM 
values of 13.01, 9.88, 7.18, and 4.87 aremin for the 70, 100, 143, 
and 217 GHz maps, respectively. 

The FWHM of the foreground cleaned products is 5 aremin 
for the maps used in this work, namely the SEVEM, SMICA, NILC, 
and COMMANDER maps. These maps are the output of four dif¬ 
ferent component-separation algorithms. While the NILC map 
is the result of an internal linear combination technique, SMICA 
uses a spectral matching approach, SEVEM a template-fitting 
method and COMMANDER a parametric, pixel based Monte Carlo 
Markov Chain technique to project out foregrounds. We refer to 
the Planck component separation papers for details in the pro¬ 
duction of these maps (Planck Collaboration IX 2015; Planck 
Collaboration X 2015). In passing, we note that the foreground 
subtraction in those maps is not perfect, and that there exist 
traces of residuals, particularly on the smallest angular scales, re¬ 
lated to radio, dust, and thermal Sunyaev-Zel’dovich [tSZ] emis¬ 
sion. 

2 Planck’s Legacy Archive: 

http://pla.esac.esa.int/pla/aio/planckProducts.html 
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For comparison purposes, we examine 9-year data from the 
Wilkinson Microwave Anisotropy Probe (WMAP) 3 . In particu¬ 
lar, we downloaded from the LAMBDA site 4 the WMAP satel¬ 
lite 9-year W-band (94 GHz) map, with an effective FWHM of 
12.4 arcmin. 

2.2. Central Galaxy Catalogue 

We define a galaxy sample in an attempt to trace the centres 
of dark matter halos. Using as a starting point the seventh data 
release of the Sloan Digital Sky Survey (SDSS/DR7 Abazajian 
et al. 2009), the Central Galaxy Catalogue (CGC) is composed 
of 262673 spectroscopic sources brighter than r = 17.7 (the r- 
band extinction-corrected Petrosian magnitude). These sources 
were extracted from the SDSS/DR7 New York University Value 
Added Galaxy Catalogue (Blanton et al. 2005) 5 , and we have 
applied the following isolation criterion: no brighter (in r-band) 
galaxies are found within 1.0 Mpc in the transverse direction 
and with a redshift difference smaller than lOOOkms . The 
Sloan photometric sample has been used to remove all possible 
non-spectroscopic sources that may violate the isolation require¬ 
ments. By using the “photometric redshift 2” catalogue (pho- 
toz2, Cunha et al. 2009), all spectroscopic sources in our isolated 
sample with potential photometric companions within 1.0 Mpc 
in projected distance, and with more than 10% probability of 
having a smaller redshift than that of the spectroscopic object, 
are dropped from the Central Galaxy Catalogue. 

Information provided in the NYU-VAGC yields estimates of 
the stellar mass content in this sample, and this enables us to 
make a direct comparison to the output of numerical simulations 
(Planck Collaboration Int. XI 2013). We refer to this paper for 
further details on the scaling between total halo mass and stellar 
mass for this galaxy sample. 

We expect most of our CGs to be the central galaxies of their 
dark matter halos (again see Planck Collaboration Int. XI2013), 
just as bright field galaxies lie at the centres of their satellite 
systems and cD galaxies lie near the centres of their clusters. 
They are normally the brightest galaxies in their system. By ap¬ 
plying the same isolation criteria to a mock galaxy catalogue 
based on the Munich semi-analytic galaxy formation model (see 
Sec. 2.3 for more details), we found that at stellar masses above 
10 11 M 0 , more than 83 % of CGs in the mock galaxy catalogue 
are truly central galaxies. For those CGs that are satellites, we 
have checked that at log(M*/M 0 ) > 11 about two-thirds are 
brighter than the true central galaxies of their halos, while the re¬ 
mainder are fainter, but are considered isolated because they are 
more than 1 Mpc (transverse direction) from their central galax¬ 
ies (60 %) or have redshifts differing by more than 1 000 km s -1 
(40%). 

A case with stricter isolation criteria has been tested, in 
which we require no brighter galaxies to be found within 
2.0 Mpc in the transverse direction and with a redshift differ¬ 
ence smaller than 2 000 km s Applying these stricter criteria 
to the same parent SDSS spectroscopic catalogue, we end up 
with a total of 110437 CGs, out of which which 58 105 galaxies 
are more massive than 10 11 M 0 . Thus about 30 % of the galaxies 
with log(M*/M Q ) >11 have been eliminated from the sample. 
This new sample of CGs has a slightly higher fraction of cen¬ 
trals, reaching about 87 % at log(M»/M 0 ) > 11. The improve¬ 
ment is small because (as we have checked) most of the satel¬ 

3 WMAP: http://map.gsfc.nasa.gov 

4 LAMBDA: http://lambda.gsfc.nasa.gov 

5 NYU-VAGC: http://sdss.physics.nyu.edu/vagc/ 


lite galaxies in the 1 Mpc sample are brighter than the central 
galaxies of their own halos, and with the stricter criteria they are 
still included. Considering a balance between the sample size, 
which directly affects the signal-to-noise ratio in our measure¬ 
ments, and the purity of central galaxies, we choose the 1 Mpc 
isolation criteria and the corresponding CGC sample in our anal¬ 
ysis from here on. This sample amounts to 262673 sources over 
the DR7 footprint (about 6 300 deg 2 , /j ky = 0.15). 

2.3. Numerical simulations 

In order to compare measurements derived from observations to 
theoretical predictions, we make use of two different numeri¬ 
cal simulations of the large-scale-structure. We first look at the 
combination of two hydrodynamical simulations, combining a 
constrained realization of the local Universe and a large cos¬ 
mological simulation covering 1200/T 1 Mpc. The simulations 
were performed with the GADGET-3 code (Springel et al. 2001; 
Springel 2005), which makes use of the entropy-conserving for¬ 
mulation of smoothed-particle hydrodynamics (SPH) (Springel 
& Hernquist 2002). These simulations include radiative cool¬ 
ing, heating by a uniform redshift-dependent UV background 
(Haardt & Madau 1996), and a treatment of star formation 
and feedback processes. The latter is based on a sub-resolution 
model for the multiphase structure of the interstellar medium 
(Springel & Hernquist 2003) with parameters that have been 
fixed to obtain a wind velocity of around 350 km s” 1 . We used the 
code SMAC (Dolag et al. 2005) to produce full sky maps from the 
simulations. For the innermost shell (up to 90 h~ l Mpc from the 
observer) we use the simulation of the local Universe, whereas 
the rest of the shells are taken from the large, cosmological box. 
We construct full sky maps of the thermal and kinetic SZ signals, 
as well as halo catalogues, by stacking these consecutive shells 
through the cosmological boxes taken at the evolution time cor¬ 
responding to their distance. In this way, the simulations reach 
out to a redshift of 0.22 and contain 13 058 objects with masses 
above 10 14 M G . Provided that our CGC lies at a median redshift 
of 0.12, this combined simulation should be able to help us in¬ 
terpret the clustering properties of peculiar velocities of highly 
biased, massive halos, on the largest scales. For more details on 
the procedure for using the simulations see Dolag & Sunyaev 
(2013). This simulated catalogue will hereafter be referred to as 
the CLUSTER catalogue. 

We have also used a mock galaxy catalogue, based on the 
semi-analytic galaxy formation simulation of Guo et al. (2013), 
which is implemented on the very large dark matter Millennium 
simulation (Springel et al. 2005). The Millennium simulation 
follows the evolution of cosmic structure within a box of side 
500 hr 1 Mpc (comoving), whose merger trees are complete for 
subhalos above a mass resolution of 1.7 X 10 10 lr l M 0 . Galaxies 
are assigned to dark matter halos following the model recipes de¬ 
scribed in Guo et al. (2011). The rescaling technique of Angulo 
& White (2010) has been adopted to convert the Millennium 
simulation, which is originally based on WMAP-1 cosmology, 
to the WMAP-7 cosmology. The galaxy formation parameters 
have been adjusted to fit several statistical observables, such as 
the luminosity, stellar mass, and correlation functions for galax¬ 
ies at z = 0. 

We project the simulation box along one axis, and assign ev¬ 
ery galaxy a redshift based on its line of sight (LOS) distance 
and peculiar velocity, i.e., parallel to the LOS axis. In this way 
we can select a sample of galaxies from the simulation using 
criteria exactly analogous to those used for our CGC based on 
SDSS. A galaxy is selected if it has no brighter companions 
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within 1 Mpc in projected distance and 1 0001<m s 1 along the 
LOS. This sample of simulated galaxies will be referred as the 
GALAXY catalogue. 

3. Methodology 

We will search for kSZ signatures in Planck temperature maps 
by implementing two different statistics on the CMB data. The 
first statistic aims to extract the kSZ pairwise momentum by fol¬ 
lowing the approach of Hand et al. (2012), which was inspired 
by Groth et al. (1989) and Juszkiewicz et al. (1998). The second 
statistic correlates the kSZ temperature anisotropies estimated 
from Planck temperature maps with estimates of the radial pecu¬ 
liar velocities. These velocity estimates are obtained after invert¬ 
ing the continuity equation relating galaxy density with peculiar 
velocities, as suggested by DeDeo et al. (2005) and Ho et al. 
(2009). The particular inversion methods applied to our data are 
described in Kitaura et al. (2012a, hereafter K12). Throughout 
this work we use the HEALPix software (Gorski et al. 1998) 6 to 
deal with the CMB maps. 

3.1. The pairwise kSZ momentum estimator 



0 200 400 600 800 1000 1200 

Distance r [A . -1 Mpc] 


Ferreira et al. (1999) developed an estimator for pairwise mo¬ 
mentum with weights depending only on line-of-sight quantities. 
With this motivation we use a similar weighting scheme on our 
CMB maps. The pairwise momentum estimator combines infor¬ 
mation on the relative spatial distance of pairs of galaxies with 
their kSZ estimates in a statistic that is sensitive to the gravita¬ 
tional infall of those objects. Specifically we use 


Pksz(r) = - 


YjKjiSTj ~ 5Tj) Ci tJ 
2i<; c y 


( 2 ) 


where the weights cy are given by 


ri + r, 


Ci.i = r L 


(n - rj)( \ + cos 6) 

2 y rf + rj - 2r,rj cos 6 


(3) 


Here r, and rj are the vectors pointing to the positions of the ith 
and /th galaxies on the celestial sphere, r, and r ; are the comov¬ 
ing distances to those objects, and ry = r, - rj is the distance 
vector for this galaxy pair. The hat symbol (r) denotes a unit 
vector in the direction of r, and 0 is the angle separating r, and 
rj. We also note that r, ; - means (r; - rj)/\r, - r,j. The sum is over 
all galaxy pairs lying a distance ry falling in the distance bin as¬ 
signed to r. The quantity ST; denotes a relative kSZ temperature 
estimate at the position of the /th galaxy: 


ST, = T AP (hi) - T AP (zi,cr z ). 


(4) 


In this equation, the symbol Tap(«/) corresponds to the kSZ am¬ 
plitude estimate obtained at the angular position of the /th galaxy 
with an aperture photometry (AP) filter, while T A p(zi,cr z ) de¬ 
notes the average kSZ estimate obtained from all galaxies after 
weighting by a Gaussian of width cr z centred on z,: 


T A p(Zi,cr z ) = 


Ti j T A p(hj) exp ( _ ^~r") 

vn / (Zj ZjP \ 

Zj ex p(— d") 


(5) 


In this AP approach one computes the average temperature 
within a given angular radius 6, and subtracts from it the av¬ 
erage temperature in a surrounding ring of inner and outer radii 


6 http://healpix.jpl.nasa.gov 


Fig. 1. Absolute value of T A p(z(r), cr z ) versus distance to the ob¬ 
server after obtaining AP kSZ temperature estimates on the real 
positions of CGs (black line) and on a rotated CG configuration 
(red line). The green line provides the theoretical expectation for 
the rms of T A p(z(r), cr z ), and the blue histogram displays the ra¬ 
dial distribution of the CGs (see text for further details). 


with / > 1. In our case we use / = V2 (see, e.g., Planck 
Collaboration Int. XIII 2014). If we knew accurately the spa¬ 
tial gas distribution on these scales there would clearly be room 
for more optimal approaches, like a matched filter technique (Li 
et al. 2014). But on these relatively small scales the density pro¬ 
file in halos seems to be significantly less regular than the pres¬ 
sure profile (see, e.g., Arnaud et al. 2010), and this adds to the 
fact that gas outside halos, not necessarily following any given 
profile, also contributes to the signal. Furthermore, we would 
like to conduct an analysis that is as model independent as pos¬ 
sible, and thus blind to any specific model for the gas spatial 
distribution. 

The correction by a redshift-averaged quantity introduced 
in the last equation above was justified by Hand et al. (2012) 
through the need to correct for possible redshift evolution of the 
tSZ signal in the sources. In practice, this correction allows us to 
measure relative changes in temperature anisotropies between 
nearby galaxy pairs after minimizing other noise sources, such 
as CMB residuals. 

We next conduct an analysis into the motivation behind this 
z-dependent monopole correction. In Fig. 1 we display the am¬ 
plitude of |7A.p(z(r),<x-)[ for cr, - 0.01 versus the comoving 
distance to the observer r. The value of the aperture chosen in 
this exercise is 8 arcmin. The black line refers to the real posi¬ 
tion of the CGs, while the red line refers to a rotated-on-the- 
sky configuration of the CGs. Since we are plotting absolute 
values of f A p{z, cr z ), the solid (dashed) parts of these lines re¬ 
fer to positive (negative) values of 7ap(z, cr,). The amplitude of 
these curves is close to the solid green line, which corresponds 
to the theoretical prediction for the rms of t A p{z, cr z ) under the 
assumption that the AP temperature estimates are dominated 
by CMB residuals. That is, the amplitude of the green curve 
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equals cr[5T}l y/Nc gM, where cr[6T] - 20yuK is the rms of the 
8 arcmin AP kSZ estimates and Ac gM is the number of CGs 
effectively falling under the redshift Gaussian window given in 
Eq. (5). While cr[5T] - 20 pK is computed from real data, its 
amplitude is very close to the theoretical predictions of Fig. 6 in 
Hernandez-Monteagudo et al. (2006b), for an angular aperture of 
8 arcmin after considering the CMB exclusively. We can see that 
the amplitude of 7’ A p(z, cr z )\ is lowest for those distances where 
the number density of CGs is highest (as displayed, in arbitrary 
units, by the blue histogram in Fig. 1). 

Therefore Fig. 1 shows not only that AP kSZ estimates are 
dominated by CMB residuals, but also that the difference in 
Tap(z, cr z ) for CG pair members lying 50-100A -1 Mpc away 
in radial distance typically amounts to few times 0.01 pK. We 
shall show below that this amplitude is not completely negli¬ 
gible when compared to the typical kSZ pairwise momentum 
amplitude between CGs at large distances, and thus subtracting 
Ta?(z, cr z ) becomes necessary. However, this should not be the 
case when cross-correlating AP kSZ measurements with esti¬ 
mates of radial peculiar velocities (see Sect. 3.2 below), since in 
this case the residual 7ap(z, <r ~) will not be correlated with those 
velocities and should not contribute to this cross-correlation. 

Like Hand et al. (2012), our choice for cr z is cr z = 0.01, al¬ 
though results are very similar if we change this by a factor of 
2. Adopting higher values introduces larger errors in the pecu¬ 
liar momentum estimates, and smaller values tend to suppress 
the power on the largest scales; we adopt cr z — 0.01 as a com¬ 
promise value. It is worth noting that any effect giving rise to a 
pair of 6Tj , 6Tj, whose difference is not correlated to the relative 
distance of the galaxies r, - r ; , should not contribute to Eq. (2). 


3.2. Recovering peculiar velocities from the CGC galaxy 
survey 


There is a long history of the use of the density field to generate 
estimates of peculiar velocities (e.g., Dekel et al. 1993; Nusser 
& Davis 1994; Fisher et al. 1995; Zaroubi et al. 1995, for some 
early studies). DeDeo et al. (2005) and Ho et al. (2009) first sug¬ 
gested inverting the galaxy density field into its peculiar velocity 
field in the context of kSZ studies. In this section, we use the ap¬ 
proach of K12 to obtain estimates of the peculiar velocity field 
from a matter density tracer. The work of K12 is based upon the 
analysis of the Millennium simulation, where full access to all 
dark matter particles was possible. In our case the analysis will 
obviously be limited to the use of the CGC catalogue, and this 
unavoidably impacts the performance of the algorithms. 

Inspired by K12, in our analysis we conduct three different 
approaches to obtain the velocity field. The first one is a simple 
inversion of the linear continuity equation of the density con¬ 
trast obtained from the galaxies in the CGC: this will be here¬ 
after called the LINEAR approach. It makes use of the continuity 
equation at linear order. 


dd(x) 

dt 


+ Vv(jc) = 0, 


( 6 ) 


where v(jc) is the peculiar velocity field and d(x) is the matter 
density contrast. In our case, however, we observe the galaxy 
density contrast 6 s (x), and thus the amplitude of the estimated 
velocity field is modulated by the bias factor b relating <5 s (jc) and 
6(x), 6 g (x) = bd{x). The bias factor b is assumed to be constant 
on the scales of interest. 

The second approach performs the same inversion of the lin¬ 
ear continuity equation, but on a linearized estimate of the den¬ 
sity field. This linearized field is obtained after computing the 


natural logarithm of unity plus the density contrast of the galaxy 
number density field and subsequently removing its spatial aver¬ 
age (Neyrinck et al. 2009), 

d LOG (jc) = In (1 + <5 g (*)) - (In (1 + <5 g )>s pa tiai. (7) 

In this expression, <5 s (jc) denotes again the density contrast of 
the galaxy number density at position x. This approach will be 
referred to as LOG-LINEAR. Finally, second-order perturbation 
theory (see, e.g., Bouchet et al. 1995) was applied on this lin¬ 
earized field, yielding a third estimate of the peculiar veloc¬ 
ity field (the LOG-2LPT approach). We refer to K12 for de¬ 
tails on the implementations of the three approaches. For the 
sake of simplicity, we discuss results for the LINEAR approach, 
and leave the corresponding discussion of the two other meth¬ 
ods for an appendix. We have also tried using the full SDSS 
spectroscopic sample (rather than the CGC) on the same vol¬ 
ume to recover the peculiar velocity field; the full galaxy sample 
should be a better proxy for the dark matter density field than the 
CGC. However, after testing this with an enlarged version of the 
GALAXY catalogue (by considering all halos above a mass thresh¬ 
old of 1O 1O ' 8 M 0 ), we obtain negligible differences with respect 
to the mock CG. Likewise, we obtain practically the same re¬ 
sults when using the real CGC and the full spectroscopic sample. 
Therefore, for the sake of simplicity, we restrict our analysis to 
the CGC. All these methods make use of FFTs requiring the use 
of a 3D spatial grid when computing galaxy number densities. 
We choose to use a grid of 128 3 cells, each cell being 4 h~ l Mpc 
on a side. This cell size is well below the scales where the typical 
velocity correlations are expected (above 40 !r x Mpc), and com¬ 
parable with the positional shifts induced by the redshift space 
distortions (about 3 A 1 Mpc for a radial velocity of 300 km s 1 ). 

These distortions will be ignored hereafter, since they affect 
scales much smaller than those of interest in our study (roughly 
20 A -1 Mpc and above). We note that these distortions can be cor¬ 
rected in an iterative fashion (see the pioneering work of Yahil 
et al. 1991 within the linear approximation; Kitaura et al. 2012c 
for the lognormal model, or Kitaura et al. 2012b including non¬ 
local tidal fields). Nevertheless, our tests performing such kinds 
of correction on the mock catalogue yield a very minor improve¬ 
ment on the scales of interest. 

After aligning the X- and Z-axes of the 3D grid with the 
zero Galactic longitude and zero Galactic co-latitude axes, re¬ 
spectively, we place our grid at a distance vector R\ vn = 
[-300, -250,150] h 1 Mpc from the observer. This position 
vector locates the corner of the 3D grid that constitutes the origin 
for labelling cells within the box. This choice of /? hl)X is moti¬ 
vated by a compromise between having as many grid cells in the 
CGC footprint as possible, and keeping a relatively high galaxy 
number density. Placing the 3D grid at a larger distance would 
allow us to have all grid cells inside the CGC footprint, but at the 
expense of probing distances where the galaxy number density 
is low (due to the galaxy radial selection function being low as 
well). Our choice for /? hox results in about 150000 CGs being 
present in the box, and about 82 % of the grid cells falling inside 
the CGC footprint. 

The three approaches will provide estimates of the peculiar 
velocity field in each grid cell, v lec (jc). From these, it is straight¬ 
forward to compute the radial component as seen by the ob¬ 
server, vj“(jc), and to assign it to all galaxies falling into that 
grid cell. In Fig. 2 arrows show the LINEAR reconstruction of 
the ^-component of the peculiar velocity field from the CGC, for 
a single z-slice of data. The coloured contour displays the galaxy 
density contrast distribution over the same spatial slice. 
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Fig. 2. Reconstruction of the ^-component of the peculiar ve¬ 
locity field (arrows) in a narrow slice of the grid containing the 
CGC over the corresponding density contrast contour 2D plot. 
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Fig. 3. Correlation coefficient of the recovered line-of-sight ve¬ 
locities with the actual ones in our GALAXY mock catalogue. 
Solid, dotted, and dashed lines refer to the LINEAR, LOG- 
LINEAR, and LOG-2LPT approaches, respectively. The red 
lines consider the ideal scenario without any sky mask or se¬ 
lection function, while the black ones are for the same sky mask 
and selection function present in the real CGC. 


The methodology outlined in K12 was conducted in the ab¬ 
sence of any sky mask or selection function. In our work we 
address these aspects of the real data by means of a Poissonian 
data augmentation approach. In a first step, all grid cells falling 
outside the CGC footprint are populated, via Poissonian realiza¬ 
tions, with the average number of galaxies dictated by the CGC 
radial selection function computed from cells inside the foot¬ 
print. In this way we fill all holes in the 3D grid. Following ex¬ 
actly the same procedure, we next radially augment the average 
number of counts in cells in such a way that the radial selection 
function of the resulting galaxy sample is constant. In this way 
we avoid radial gradients that could introduce spurious velocities 
along the line of sight. 

In order to test this methodology, we use the GALAXY mock 
catalogue, in two different scenarios. The first scenario is an 
ideal one, with no mask or selection function: we apply the three 
approaches on a box populated with our mock catalogue, and 
compare the recovered radial velocities with the real ones pro¬ 
vided by the catalogue. When performing this comparison, we 
evaluate a 3D grid for the original galaxy radial velocity v °“ s by 
assigning to each cell the radial velocity of the galaxy falling 
nearest to the cell centre, that is, we adopt a nearest particle 
(NP) method. For different radial bins of the k wavevector, we 
evaluate the correlation coefficient (r>) computed from the cross¬ 
power spectrum between the recovered and the original radial 
velocity fields and their respective auto-spectra: 

p orig,rec(£) 

n = . - . (8) 

yjP on S’ °ng (k) P rec ’ rec (A) 

In this expression, P x - Y (k ) stands for (v(£)^ s (v(£)^ s )*), that is, 
the radially averaged power spectrum of the line of sight velocity 
modes. The superscripts { X , Y ) = jorig,rec) stand for “original” 
and “recovered” components, respectively, and the asterisk de¬ 
notes the complex conjugate operation. 

In Fig. 3 the red lines display the correlation coefficients for 
the ideal scenario. Solid, dotted, and dashed lines refer to the 
LINEAR, LOG-LINEAR, and LOG-2LPT approaches, respec¬ 
tively. On large scales (low k modes), the three approaches pro¬ 
vide correlation coefficients that are very close to unity, while 
they seem to lose information on small scales in the same way. 
We note that a direct comparison to the results of K12 is not 
possible, since, in our case, we do not use the full dark matter 
particle catalogue, but only a central halo one. 

After the real mask and selection function obtained from 
the CGC is applied to the GALAXY mock catalogue, then the 
recovery of the peculiar velocities from the three adopted ap¬ 
proaches worsens considerably. The black lines in Fig. 3 dis¬ 
play, in general, much lower correlation levels than the red ones. 
The LINEAR approach seems to be the one that retains most 
information on the largest scales, and it out-performs the LOG- 
LINEAR and the LOG-2LPT methods on all scales. We hence 
expect the LINEAR approach to be more sensitive to the kSZ 
effect than the other two methods, particularly on the largest 
scales. 

Once the velocity inversion from the CGC density contrast 
has been performed, we compute the spatial correlation func¬ 
tion between the recovered velocities and the kSZ temperature 
anisotropies, 

w r,v (r) - (STivgQcjiXjir), (9) 

where 6Tj is estimated as in Eq. (4), and the ensemble average 
is obtained after running through all galaxy pairs {/, /'} lying a 
distance r away. 
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3.3. Template fitting 

When studying the measurements of the kSZ pairwise momen¬ 
tum and kSZ momentum-v^ correlation, we perform fits to es¬ 
timates obtained from our numerical simulations. That is, we 
minimize the quantity 

X~ = 

Uj 

where the indexes i, j run over different radial bins, w x (r,) is 
the measured quantity in the zth radial bin (with X either denot¬ 
ing kSZ pairwise momentum or the kSZ temperature-recovered 
velocity correlation), and w x Slm (r,) refers to its counterpart mea¬ 
sured in the numerical simulation. The symbol C, 7 denotes the 
i j component of the covariance matrix C that is computed from 
the Planck maps after estimating w x for null positions where 
no SZ effect is expected 7 * * . This minimization procedure provides 
formal estimates for the amplitude A x and its associated errors: 

y C-'w X ’" m ( rj ) 

A x — J _ bl _( 11 ) 

w x ’ sim (r,) Crj W X ' sim (r 7 ) ’ 


2/j w* sim (r,) C~j w x > ™( rj ) ’ 

Most of the information is located at short and intermediate 
distances, where the estimated statistic w x (r) differs most from 
zero, as can be seen in Figs. 4 and 5. We also test the null hy¬ 
pothesis, that is, we measure the x 2 statistic (defined in Eq. 10 
above) for the particular case of A x = 0 and estimate the sig¬ 
nificance level at which such a value (denoted by ^ u]] ) is com¬ 
patible with this null hypothesis. In these cases, we quote the 
significance as the number of <x with which the null hypothesis 
is ruled out under Gaussian statistics. The array of distance bins 
on which the covariance matrix in Eq. (10) is computed is cho¬ 
sen evenly in the range 0—150/z 1 Mpc. However, we consider 
only three separate points centred upon 16, 38, and 81 h~ 1 Mpc 
when computing (conservative) statistical significances. In this 
way we minimize correlation among radial bins to ensure that 
the inversion of the covariance matrix is stable. We have checked 
that, for the adopted set of distance bins, random variations at the 
level of 10% of the measured w x (r ) do not compromise the sta¬ 
bility of the recovered significance estimates. In other words, we 
explicitly check that 10 % fluctuations on the measured vi^(r) in¬ 
troduces fluctuations at a similar level in the x 2 estimates (more 
dramatic changes in the x 2 estimates would point to singular or 
quasi-singular inverse covariance matrices). 

4. Results 

4. 1. The kSZ pairwise momentum 

As mentioned above, the CGC was used previously in Planck 
Collaboration Int. XI (2013) to trace the tSZ effect versus stellar 
mass down to halos of size about twice that of the Milky Way. 
In this case we use the full CGC to trace the presence of the kSZ 
signal in Planck data, since our attempts with the most massive 
sub-samples of the CGC yield no kSZ signatures. While the tSZ 
effect is mostly generated in collapsed structures (Hernandez- 
Monteagudo et al. 2006a) because it traces gas pressure, the kSZ 

7 These null positions on the Planck maps correspond to rotated or 

displaced positions with respect to the original location of the CGs, as 

will be explained in the next section. 


effect instead is sensitive to all baryons, regardless of whether 
they belong to a collapsed gas cloud or not. Thus it is expected 
that not only will the virialized gas in halos contribute to the kSZ 
signal, but also the baryons surrounding those halos and moving 
in the same bulk flow. Rather than using a particular gas density 
profile, we choose not to make any assumption about the spatial 
distribution of gas around CGs. 

We thus adopt an AP filter of varying apertures around the 
positions of CGC galaxies. The minimum aperture we consider 
is close to the resolution of Planck (5 arcmin), and we search 
for signals using increasing apertures of radii 5, 8, 12, and 
18 arcmin. This is equivalent to probing spheres of radius rang¬ 
ing from 0.5 up to 1.8 /z _1 Mpc (in physical units) around the 
catalogue objects. The result of the kSZ pairwise momentum es¬ 
timation on the raw 217 GHz HFI map is displayed in the top 
row of Fig. 4 for the four apertures under consideration. The 
recovered momenta for the raw frequency maps (displayed by 
black circles) provide some evidence for kSZ signal for aper¬ 
tures smaller than 12 arcmin: below a distance of 60 hr 1 Mpc 
all points are systematically below zero, some beyond the 2<x 
level. Although the test does not yield significant values 
(at 0.4cr,0.3cr,0.7cr and -1.2cr for 5, 8, 12, and 18arcmin 
apertures, respectively), the fits to the GALAXY peculiar momen¬ 
tum template (displayed by the green solid line in the second- 
from-the-left panel in the top row) yield S/N = A Pksz /cr Pksz — 
1.7,1.4,1.9 and 0.1 for 5, 8, 12, and 18arcmin apertures, re¬ 
spectively. We obtain similar levels for the peculiar momen¬ 
tum template obtained from the CLUSTER simulation (red line in 
second-from-the-left panel in top row). When using the SEVEM 
clean map, these significances go beyond the 2 cr level: S/N 
= A Pksz /cr Pksz = 2.5,1.8,2.2 and 0.0 for 5, 8, 12, and 18 arcmin 
apertures for the GALAXY-derived template, respectively. If we 
repeat the analysis at 12 arcmin aperture for the SMICA NILC 
and COMMANDER maps, we obtain S/N= A Pksz /cr Pksz — 2.1,2.2, 
and 2.1, respectively. 

We note at this point that the behaviour of the kSZ evidence 
for the CGs as displayed in Fig. 4 is significantly different to 
what is found in Fig. 5 for the CLUSTER simulation with halos 
in the mass range (1-2) X 10 14 h 1 M Q . In this simulation, most 
of the kSZ signal is coming from the halos themselves, and thus 
increasing the aperture to values larger than the virial radius of 
the clusters results in a dilution of the kSZ pairwise momentum 
amplitude. For real data, we find that the amplitude of the signal 
does not show significant changes when increasing the aperture 
from 5 arcmin up to 12 arcmin. It is worth adding now that the 
black filled triangles in Fig. 5 represent the peculiar momentum 
measured on the GALAXY mock catalogue after assigning to these 
galaxies a tSZ temperature fluctuation following the tSZ versus 
mass scaling measured in the CGC (Planck Collaboration Int. XI 
2013). This tSZ amplitude assumes an observing frequency of 
fobs = 100 GHz, and hence provides a conservative estimate of 
the tSZ contamination, which turns out to be at the level of -3 x 
10 -3 pK, i.e., about a factor of 30 below the measured amplitude 
on the real CG sample. 

We next investigate the spectral stability of the signal for 
12 arcmin apertures. This choice is motivated by the fact that the 
LFI 70 GHz and the HFI 100 and 143 GHz channels have lower 
angular resolution, comparable or bigger than 8 arcmin, and this 
may compromise the comparison with the 217 GHz channel. 
This analysis is shown in the bottom row of Fig. 4. The signal 
obtained for the raw 217 GHz channel is found to be remarkably 
similar to what we obtain at 143, 100 GHz in HFI: the fits to the 
kSZ pairwise momentum template from the GALAXY simulation 
gives S/N = A Pksz /cr Pksz = 2.2 and 2.1, respectively. The LFI 
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Fig. 4. Computation of the kSZ pairwise momentum for the CGC sample. The top row, from left to right, displays the results for 
different aperture choices on the raw HFI 217 GHz map, namely 5, 8, 12, and 18 arcmin. The top row also shows the analysis for 
the foreground-cleaned SEVEM map, displayed with blue squares. The fit to the pairwise momentum templates from the CLUSTER 
catalogue is also displayed by the solid red line, and the GALAXY catalogue by the green line. The bottom row presents the results at 
a fixed aperture of 12 arcmin for different frequency maps, including WMAP-9 W-band data (red squares). 


70 GHz channel, on the other hand, seems to show a much flatter 
pattern, and this could be due to a larger impact of instrumental 
noise: the tSZ should also be present at 100 and 143 GHz, and, 
as expected, gives rise to negligible changes. The effective full 
width half maximum (FWHM) for the 70 GHz channel is close 
to 13 arcmin, and this may also contribute to the difference with 
respect to the other HFI channels. 


In the bottom row, second panel of Fig. 4, the red squares 
provide the measurement obtained from the cleaned W-band 
map of WMAP-9. Surprisingly, the level of anti-correlation for 
distances below 30 h~ l Mpc appears higher for WMAP data than 
for the Planck channels: the first three radial bins lie at the 2.3- 
3.3 cr level, and a fit to the GALAXY template of the pairwise 
momentum yields S/N = A Pksz /cr Pksz = 3.3 as well. The angu¬ 
lar resolution of the W-band in WMAP is close to that of the 
LFI 70 GHz channel, and the non-relativistic tSZ changes by 
less than 13 % between those channels. If any other frequency- 
dependent contaminants (which should be absent in the WMAP 
W-band) are responsible for this offset, then they should also 
introduce more changes at the other frequencies. Thus the rea¬ 
son for the deeper anti-correlation pattern found in the W-band 
of WMAP-9 remains unclear. Overall, we conclude that the 
Planckresults for the pairwise peculiar momentum are compat¬ 
ible with a kSZ signal, at a level ranging between 2 and 2.5 cr. 
The output from WMAP-9 W band is also compatible with a 
kSZ signal, with a statistical significance of 3.3 cr. 


4.2. Cross-correlation analysis with estimated peculiar 
velocities 

After reconstructing the CGC density field into a peculiar ve¬ 
locity field in a 3D grid, as explained in Sect. 3.2, we compute 
the cross-correlation function w T ’ l 4«(r) for the three inversion ap¬ 
proaches using the Planck HFI 217 GHz band, and the SEVEM , 
SMICA , NILC , and COMMANDER maps (see Fig. 6). Here we are 
using kSZ temperature estimates obtained with an 8 arcmin aper¬ 
ture. The three inversion approaches (namely, LINEAR, LOG- 
LINEAR, and LOG-2LPT) may suffer from different biases that 
impact the amplitude of the recovered velocities. Thus, in order 
to avoid issues with the normalization of the recovered veloc¬ 
ities and to ease the comparison between different approaches, 
we normalize the recovered velocities v[“-s by their rms before 
computing the cross-correlation function. In any case we must 
keep in mind the different responses of the three approaches on 
different k-scales (see Fig. 3), thus giving rise to different corre¬ 
lation structures. For the sake of clarity, in this section we show 
results for the LINEAR method only, while in the Appendix we 
discuss the results for the other two reconstruction methods. 

In Fig. 6, filled coloured circles display the results obtained 
after using the 6T kSZ temperature estimates obtained at the real 
positions of the CGs on the CMB maps. The null tests were 
obtained by computing the cross-correlation of the normalized 
recovered velocities (v[“) with temperature estimates (6T) ob¬ 
tained on 50 rotated 8 positions on the CMB maps. In those 
cases, the 6T s were computed for positions obtained after ro- 


8 For each rotation/displacement we use a step of three times the 
aperture radius adopted. 
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Fig. 5. (Filled circles): Measured kSZ pairwise momentum for 
the kSZ map derived from the CLUSTER simulation after consid¬ 
ering different radial apertures on a subset of clusters in the mass 
range (1-2) xlO 14 h~ ] M e . (Black triangles): Pairwise momen¬ 
tum computed from the GALAXY mock catalogue (filled triangles) 
after assigning these galaxies a tSZ amplitude following a mass 
scaling inspired by the tSZ measurements of the CGC given in 
Planck Collaboration Int. XI (2013). 


tating the real CG angular positions in Galactic longitude. The 
results from each of these null rotations are displayed by dotted 
lines in Fig. 6, and their average is given by the thick, dashed 
line (which lies close to zero at all radii). Both the error bars and 
the covariance matrix of the correlation function were obtained 
from these null realizations. 

Results at zero-lag rotation (i.e., the real sky) lie far from 
the distribution of the null rotations. There are several individual 
distance bins lying more than 3 cr (up to 3.8 cr on the 16 h~ l Mpc 
distance bin for the raw 217 GHz frequency map), and the X^uU 
tests rule out the null hypothesis typically at the level of 2.1- 
2.6 cr for the clean maps, and at 3.2 cr for the raw 217 GHz 
frequency map. Likewise, when fitting the observed correlation 
function to the w l,v i<» (r) correlation function obtained from the 
GALAXY simulation, we obtain values of S/N = A H ' T '' /cr"' 1 '' = 
3.0-3.2 for the clean maps, while for the HFI 217 GHz raw map 

Tv Tv 

we find S/N = A w ' /cr w ' = 3.8. We explain this apparent mis¬ 
match below. A clear, large-scale correlation pattern, extending 
up to about 80 h~ l Mpc, is found in the data. A complemen¬ 
tary systematic test can be conducted by computing the cross¬ 
correlation function of the kSZ temperature fluctuations ( 6T ) 
with shuffled estimates of the recovered line-of-sight peculiar ve¬ 
locities, i.e., to each CG we assign a v[ ec estimate corresponding 
to a different, randomly selected CG. The result of performing 
this test for the LINEAR approach is displayed in Fig. 7, and 
shows that the correlation found between the 5T s and the recov¬ 
ered velocities clearly vanishes for all shuffled configurations. 
By shuffling the recovered velocities we are destroying their co¬ 
herent, large-scale pattern, which couples with large angle CMB 
residuals, and generates most of the uncertainty in the measured 



Fig. 6. Measured cross-correlation function between the kSZ 
temperature estimates and the recovered radial peculiar veloc¬ 
ities, w T,v '"'(r), according to the LINEAR approach. The re¬ 
covered velocities are divided with their rms dispersion cr v = 
310 km s _1 , so kSZ temperature estimates are correlated to a 
quantity of variance unity. This plot corresponds to an aper¬ 
ture of 8 arcmin. Filled coloured circles correspond to iv T '' : L(r) 
estimates from different CMB maps (SEVEM, SMICA, NILC, 
COMMANDER, and the HFI 217 GHz map). The dotted lines dis¬ 
play the null estimates obtained after computing kSZ tempera¬ 
ture estimates for rotated positions on the SEVEM map, and the 
thick dashed line displays the average of the dotted lines. Error 
bars are computed from these null estimates of the correlation 
function. The solid line provides the best fit of the SEVEM data to 
the theoretical prediction for w l A '<»(r) obtained from the GALAXY 
mock catalogue. These predictions are obtained using only a rel¬ 
atively small number of mock halos, and hence their uncertainty 
must be considered when comparing to the data. 


cross-correlation. This explains the smaller error bars in Fig. 7 
when compared to Fig. 6. 

We further study the dependence of the measured w 1 ,V los W 
correlation function on the aperture radius and show the results 
in Fig. 8. We again restrict ourselves to the LINEAR approach 
and the SEVEM foreground-cleaned map, and error bars are com¬ 
puted from the null rotations (rather than from shuffling esti¬ 
mates of v™). We display results for four apertures ranging from 
5 to 18 arcmin: the lower end of this range is given by the angu¬ 
lar resolution of the map, while for the higher end the kSZ signal 
is found to vanish. We find that, as for the kSZ peculiar momen¬ 
tum, the kSZ amplitude at 8 arcmin aperture is very close to the 
one found at 5 arcmin, and it is still significant at 12 arcmin aper¬ 
ture. We then vary the aperture from 5 to 26 arcmin and calculate 
the corresponding x 2 values, as shown in Fig. 9. These results 
are for the SEVEM foreground-cleaned map, and provide another 
view of the angular extent of the signal. For the kSZ peculiar 
momentum, we have shown in Fig. 4 that there is kSZ evidence 
for apertures as large as 12 arcmin, in good agreement with what 
we find now for the w T ’'’I™(r) correlation function. While this 
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Fig. 7. Same as in Fig. 6, but with the null correlation functions 
(displayed by dotted lines) being estimated after assigning the 
recovered linear velocity estimate of a given CG to any other 
CG in the sample (that is, by “shuffling”). Only the LINEAR 
approach is displayed in this plot. 


statistic seems to have higher significance than the pairwise mo¬ 
mentum, for both statistics we find consistently that most of its 
signal is again coming from gas not locked in the central regions 
of halos, but in the intergalactic medium surrounding the CG 
host halos. 

We also test the consistency of our results with respect to fre¬ 
quency on raw Planck raw maps. For that we use the HFI chan¬ 
nels ranging from 100 GHz up to 353 GHz and now we fix the 
aperture at 8 arcmin. Since the 100 and 143 GHz frequency maps 
have angular resolution comparable or worse than 8 arcmin, we 
choose to deconvolve all HFI maps under consideration by their 
respective (approximate) Gaussian beams, and convolve them 
again with a Gaussian beam of FWHM= 5 arcmin. While this 
approach may challenge the noise levels for the HFI maps with 
coarser beams, we find that this is balanced by the large number 
of CGs on which we compute the w' -'io.fr) correlation function. 
The results of this analysis are given in Fig. 10, and show how 
the w'-'iEfr) correlation functions from 100GHz up to 217 GHz 
agree closely with each other, with no significance dependence 
on frequency (as is expected for the kSZ effect). However, the re¬ 
sult for the 353 GHz channel is pointing to a significantly higher 
amplitude of the correlation function, even if the error bars as¬ 
sociated with this map are typically 50-70% larger than for the 
100 GHz map. At the map level, the main difference between 
the 353 GHz and lower frequency channels is the significantly 
larger amount of dust and/or CIB emission in the former map. 
We have checked that, throughout the rotated configurations, the 
HFI 353 GHz map has on average no correlation with the esti¬ 
mated radial velocities of the CGs. Therefore, the excess found 
in the measured amplitude of the correlation function 

must be due to fortuitous alignment between the estimated ra¬ 
dial velocities of the CGs and the dust emission at the position 
of the CGs. This would also explain the higher amplitude of the 


Fig. 8. Dependence on the aperture radius of the (r) corre¬ 

lation function obtained in the LINEAR approach for the SEVEM 
map. Error bars are estimated from the rms of the null correlation 
functions computed from rotated kSZ temperature estimates. 


w T v 1S(r) correlation function found for the raw HFI 217GHz 
map with respect to all other clean maps, as shown in Fig. 6: the 
raw HFI 217 GHz map should still contain some non-negligible 
dust contamination when compared to the foreground-cleaned 
maps SEVEM, SMICA, NILC, and COMMANDER. 

Another potential issue is related to the fact that, as shown 
in Planck Collaboration Int. XI (2013), about 12% of the CGC 
members with stellar masses above 10 10 M Q are actually not cen¬ 
tral galaxies, but show some offset with respect to the centres of 
the host halos. This offset should cause a low bias in the kSZ 
amplitude estimation. We have simulated the impact of this off¬ 
set by assuming that 12 % of the CGs have uniform probability 
of lying within a (projected) distance range of 0 to 1 Mpc from 
the halo centre. After adopting a NFW profile for the gas distri¬ 
bution, we have found that the LOS-projected kSZ signal should 
be low biased a 10, 8, and 6% for $ap = 5,10 and 15 arcmin, 
respectively. Given the limited precision of our measurements, 
these biases are relatively small and will be ignored hereafter. 

To compare our detection with simulations, the solid line 
in Fig. 6 shows the best fit to the prediction inferred from the 
GALAXY mock catalogue. This prediction is obtained by applying 
the LINEAR velocity recovery algorithm on our GALAXY cata¬ 
logue after imposing the sky mask and the selection function 
present in the real CGC. Three pseudo-independent estimates of 
the correlation function of the recovered LOS velocity v rec and 
the real LOS velocity vi os ((viosHecXf) = w Vr “’ l ’(r)) are obtained 
after rotating the 3D grid hosting the mock GALAXY CGC, so 
that the side facing the observer is different in each case. The 
solid line corresponds to the average of these three estimates of 
w' ?i “' v (r), and the ratio to the observed correlation function can be 
interpreted as an “effective” optical depth to Thomson scattering: 
w T ’'k>s(r) - —r T w Vr “ ,v (r). We obtain tt = (1.39 + 0.46) x 10~ 4 
(i.e., at the 3 cr level) for the SEVEM map, with very similar values 
for all other foreground-cleaned maps. We defer the physical in- 
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Fig. 9. Evaluations of the^ 2 statistic of the w T ’'!“(>) correlation 
function (with respect to the null hypothesis) for different an¬ 
gular apertures. Only results for the foreground-cleaned SEVEM 
map are shown. There is evidence for kSZ signal in a wide range 
of apertures above the FWHM, with a local maximum close to 
8 arcmin. This roughly corresponds to a radius of 0.8 /z 1 Mpc 
from the CG positions at the median redshift of the sample. 

terpretation of the kSZ measurements of this paper to an external 
publication, (Hernandez-Monteagudo et al. 2015). 

5. Discussion and Conclusion 

The roughly 2.2 cr detection (varying slightly for the different 
maps) of the pairwise momentum indicates that the baryonic gas 
is comoving with the underlying matter flows, even though it 
may lie outside the virial radius of the halos. The aperture of 
8 arcmin on the CGs (placed at a median redshift of z - 0.12) 
corresponds to a physical radius of around 1 Mpc. As we show 
next, this is considerably higher than the typical virial radius of 
the CG host halos. Following the same approach as in Planck 
Collaboration Int. XI (2013), we compute the /Goo radius con¬ 
taining an average matter density equal to 200 times the critical 
density at the halo’s redshift. In Fig. 11 we display the histogram 
of the angle subtended by the /Goo values of the 150000 CGs 
placed in the 3D grid that we use for the velocity recovery. The 
red vertical solid line indicates the 8 arcmin aperture, well above 
the typical angular size of the CGC sources. 

The behaviour displayed by the measured kSZ peculiar mo¬ 
mentum in the top row panels of Fig. 4 differs significantly from 
the pattern found in Fig. 5. The fact that the measured kSZ pair¬ 
wise momentum shows a roughly constant amplitude out to an 
aperture of 12 arcmin, well above the CG virial size, signals the 
presence of unbound gas that is contributing to the measurement. 
The opposite situation is seen in Fig. 5: in this scenario most of 
the signal comes from the halos rather than from a surrounding 
gas cloud. This plot displays the kSZ peculiar momentum from a 
subset of our CLUSTER catalogue, with sources in the range (1- 
2)xl0 14 /; 1 M c .,, after considering different aperture radii. As 


FWHM eff =5 arcmin, AP = 8 arcmin 



Fig. 10. Variation of the recovered w'^'P (r) correlation function 
for the four lowest frequency raw HFI maps after an effective 
convolution by a Gaussian beam of FWHM = 5 arcmin. Error 
bars here are computed in the same way as in Figs. 6 and 8. 



9 200 [arcmin] 

Fig. 11. Histogram of the angle subtended by /Goo of the CGC 
members. 


long as the halo remains unresolved (as is the case for these sim¬ 
ulated clusters), then as we increase the aperture size the kSZ 
signal coming from the halo becomes more diluted and hence 
the amplitude decreases, contrary to what is found in Fig. 4. 

We find a similar situation in the kSZ temperature- 
peculiar velocity cross-correlation. By cross-correlating the re¬ 
constructed peculiar velocity field in a 3D box with the kSZ tem- 
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perature anisotropies, we find a 3.0 cr detection between the two 
fields for the SEVEM map, at an aperture of 8 arcmin. This again, 
corresponds to gas clouds with radius roughly 1 Mpc, about 
twice the mean /Goo radius (we find that (/Goo)cgc — 0.4 Mpc). 
Since the peculiar velocity is directly related to the underlying 
matter distribution, our result suggests that gas inside and out¬ 
side CG host halos are comoving with the matter flows. 

One way of quantifying the amplitude of our signal is to 
ask by what factor we need to scale the model-based w l ’ rec,v (r) 
in order to match our measured w T '^(r). Interpreting this scal¬ 
ing as an “effective” optical depth to Thomson scattering we find 
Tt = (1.4 + 0.5) x 1CL 4 , which is a factor of 3 larger than that ex¬ 
pected for the gas in typical CG host halos alone. This provides 
another piece of evidence that the kSZ signal found in Planck 
data is generated by gas beyond the virialized regions around 
the CGs, as opposed to the tSZ effect, which is mostly gener¬ 
ated inside collapsed structures (Hernandez-Monteagudo et al. 
2006a; Van Waerbeke et al. 2014; Ma et al. 2014) . 
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Appendix A: The three velocity reconstruction 
methods 

In this appendix we briefly compare the velocity recovery by 
means of the three approaches explained in this work, namely the 
LINEAR, the LOG-LINEAR, and the LOG-2LPT approaches. 
First of all, as is shown in Fig. 3, the LINEAR method outper¬ 
forms our implementation of the other two methods in the pres¬ 
ence of a mask and selection function (even if all three provide 
virtually identical correlation coefficients in the absence of a 
mask and selection functions). This is shown clearly in Fig. A.l. 

We apply our algorithms to invert density into velocities after 
using the three different sides of the cubic grid where we place 
our GALAXY mock catalogue. For each box side we estimate a 
recovered velocity field and compute its correlation to the real, 
underlying velocity field of the mock galaxies, w' r,!C ’ v (r). Three 
different estimations are too few to an average correlation func¬ 
tion, but we still use this average function (out of these three 
box-side estimates) when conducting the fit to real data. Thus 
the solid lines in Fig. A.l show the fits to the result of averaging 
the correlation functions w T,v through the three box orientations. 
We must keep in mind that we expect a high level of uncertainty 
in these predictions. 

These average w T l ’(r) functions show a smooth behaviour 
for the LINEAR and LOG-2LPT implementations, although 
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Fig.A.l. Measured cross-correlation function between the kSZ temperature estimates and the recovered radial peculiar velocities, 
w T l '[os(r), according to the three adopted approaches, LINEAR (left panel), LOG-LINEAR (middle panel), and LOG-2LPT (right 
panel). The velocities are normalized with their rms dispersion cr v = 230, 310, and 200 km s~ ! from left to right, respectively. We use 
an aperture of 8 arcmin when estimating kSZ temperature fluctuations. Filled coloured circles correspond to w T,v i‘*(r) estimates from 
different CMB maps (SEVEM, SMICA, NILC, COMMANDER, and the HFI 217 GHz map). The dotted lines display the null estimates 
obtained after computing kSZ temperature estimates for rotated positions on the SEVEM map, and the thick dashed line displays the 
average of the dotted lines. Error bars are computed from these null estimates of the correlation function. The solid line provides 
the best fit to the data of the theoretical prediction for w r,v '“ (r) obtained from the GALAXY mock catalogue. These predictions are 
obtained using only a relatively small number of mock halos, and hence their uncertainty must be considered when comparing to 
data. 


the LOG-LINEAR case displays a slight flattening on 50- 
80/T 1 2 3 4 5 6 7 8 9 10 Mpc scales, which is absent in the other two cases and 
is not well fit by the data. We expect correlation functions ob¬ 
tained from the GALAXY mock catalogue to be noisier than the 
real ones obtained from the CGC, since typically only about 
80000 mock galaxies survive the mask and the selection func¬ 
tion, as opposed to the 150000 surviving galaxies in the real CG 
sample. The fits of the real data to the solid lines in Fig. A.l 
yield S/N = = 3.0, 2.0, and 2.1 for the LINEAR, 

LOG-LINEAR, and LOG-2LPT approaches, respectively. 
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